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A STUDY OF MULTIPLE-T.HAK! R MODAL SURVEY TESTING* 

I . Summary 

This research has been concerned with methods of structural dynamic 
testing. The principal objective has been to examine and to assess the 
practical value of a method of multiple-shaker sinusoidal modal vibration 
testing known as Asher's method. Numerical studies which simulate the 
application of Asher's method and a unique experimental Implementation of 
the method have been completed. Another objective of the research has 
been to develop and to demonstrate with numerical simulation a quantita- 
tive method for determining from transfer function data the number of dom- 
inant modes of vibration in sinusoidal structural response. 

II. Studies of Asher's Method 

The first report on numerical simulation of Asher's method was the 
Master's thesis by Stafford [1]. Simulated modal testing on a relatively 
simple mathematical structural model was examined in considerable detail. 
Stafford's model is a cantilevered plane grid structure having five degrees 
of freedom, a pair of closely spaced modes, and hysteretic damping which 
does not couple the undamped normal modes. The next report on numerical 
simulation was the Master's thee i s by Shostak [2]. Shostak's mathematical 
models are similar Ic Stafford's but have viscous damping which does couple 
tiie undamped normal modes. The final and .lost significant report on numeri 
cal simulation was the paper by Hallauer and Stafford, which was presented 

♦The NASA Technical Officer for this grant is Mr. Robert Mlserentino, 

NASA Langley Research Center. 
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at a technical conference [3]. and published in a journal [4], The ab- 
stract for this paper follows. 

The method proposed by Asher for structural dynamic modal test- 
imj by multiple-shaker sinusoidal excitation is reviewed, and 
its theory arid application are discussed In detail. Numerical 
results from simulated modal testing on mathematical structural 
models are presented to Illustrate the strengths and weaknesses 
of the method. The characteristics of these models include 
damping which couples the normal modes and closely spaced modes. 
Numerical techniques required for implementation of the method 
are described. A procedure is suggested for replacing actual 
mechanical tuning with calculations employing transfer function 
data. 

The experimental study of Asher's method was conducted in the Struc- 
tural Dynamics Research Laboratory cf NASA Langley Research Center. The 
objective was to develop and test the software and procedures for appli- 
cation of the method with the use of the SDRL's programmable Hewlett- 
Packard 54b 1 B Fourier Analyzer System with Modal Analysis Option. The most 
detailed report on this experimental study is the Master's thesis being 
prepared by R. R. Gold, which will be submitted to VPI & SU in late 1978 
or early 1979. Also, a technical paper by R. R. Gold and W. L. Ilallauer, 

Jr. is being submitted for presentation at the 25th International Instru- 
mentation Symposium (May 7 - 10, 1979 in Anaheim, California) and publica- 
tion in the symposium proceedings. This paper is entitled "Implementation 
of Asher's Method of Modal Testing on a Fourier-Analysis/Modal-Test System", 
and its abstract follows. 

A software package for application of Asher's method has 
been Implemented on a modern Fourier analyzer system. This com- 
bination is unusual since the original form of Asher's method 
involves mul tiple-shaker, sinusoidal-dwel 1 mechanical tuning, 
whereas the more recent method of analysis-based modal testing 
employs single-point transient or random excitation, with modal 
separation attempted by curve fitting of FFT-calculated transfer 
functions. The technique described uses FFT and curve-fitting 
capabilities to produce an analytical form of multiple-shaker 


tuning. For evaluation of this technique, a model having two 
modes with almost identical natural frequencies was designed 
mathematical ly and fabricated. Evaluation with numerically 
simulated data for this model demonstrated that the technique 
works extremely well in principle. However, actual testing 
of the model produced such poor transfer function data that 
practical evaluation was not possible. 

The principal application of multiple-shaker modal testing Is in 

the difficult task of separating and correctly characterizing closely 

spaced modes of vibration. In order to assess the practical value of 

Asher's method, it has therefore been necessary to have available lelativel.v 

simple structural models with closely spaced modes. It is very difficult 

to create such models by intuition or trial and error, so a quantitative 

design method has been developed and is described in detail in the paper 

by Hallauer et al [5]. The abstract for this paper follows. 

Vehicle structures often have closely spaced moles of vibration 
within the frequency spectrum of applied loads. Such modes are 
important since they tend to be the source of vibration problems. 

In order to study the effects that closely spaced modes have on 
structural response in such situations as modal testing and self- 
cxcited vibrations, it is useful to be able to design a mathematical 
structural model having closely spaced modes. In this paper, a 
method for designing such a model is presented and illustrated 
with examples. Given a reference model with specified geometry 
and degrees of freedom, the lumped inertias and stiffnesses of 
the model are perturbed in such a manner as to force together two 
of its natural frequencies. With a slight alteration, the method 
is also applicable to the inverse problem of separating undesirable 
closely spaced modes which appear in a structural design. 

An interim report of this research (Semi-Annual Progress Report, May- 

November 1977, Research Grant NSG 1276) stated that a NASA Contractor 

Report would be published with the title "A Method for Modal Tuning by 

Multiple-Shaker Sinusoidal Excitation: Theory and Numerical Simulation". 

The CR has not been and wi 1 1 not be published. All the material that was 

to have been included in the CR is available in references [4] and Lb]. 
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111. Determining the Number of Dominant Vibration Modes 

Very significant progress has been made In the development and demon 
stration of a quantitative method lor determining from transfer function 
data the number of dominant modes. The first report on this subject 
was the Master's thesis by Franck (.6], The research was also presented 
at a technical conference in a paper by Hallauer and Franck [7]. Finally, 
a detailed and comprehensive paper by W. L. Hallauer, Jr. and A. Franck 
has been submitted for publication In late 1979 In the 48th Shock and 
Vibration Bulletin. In order to make this paper available prior to its 
publication, it has been Included as the Appendix to this report. 
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ON DETERMINING THE NUrlflER OF DOMINANT MODES 


IN SINUSOIDAL STRUCTURAL RESPONSE 
W. L. Hallauer, Jr. and A. Franck 
ABSTRACT 

This paper addresses the problem of using structural dynamic transfer 
function data to determine the number of vibration modes dominant in re- 
sponse at a given frequency. If two or more modes are closely spaced or 
if response is influenced strongly by distant modes, then the number of 
dominant modes may not be evident from examination of transfer function 
plots, and quantitative methods may be required. Two relatively simple 
methods which have been used previously are reviewed, and a more effective 
new method, called the vector-fit method, is described In detail. Applica 
tions of these methods are given with the use of numerically simulated 
transfer functions data. 


r\V<Vo,nx> WOT , ; ILto£o 


I. INTRODUCTION 


At any particular frequency of excitation, the steady-state sinusoidal 
response of a structure Is dominated by only a few of Its Indefinitely large 
number of vibration modes. The objective of nodal testing Is to measure 
specific parameters of the dominant modes such as natural frequencies, damp- 
ing values, and mode shapes. When applicable, the best method for determin- 
ing the number of dominant modes in a frequency band is simply to count 
resonance peaks on transfer function plots. In such a situation, modal param- 
eters can then be calculated rather easily by modern curve-fitting algo- 
rithms, most of which require the number of dominant modes as an input value. 
However, if two or more modes are closely spaced, or if response is influenced 
strongly by modes whose resonances are outside the freouency band of interest, 
then peak counting may not reveal the true number of dominant modes, and 
subsequent curve fitting of transfer function data may produce Incorrect 
modal parameters and/or miss modes entirely. But a quantitative method for 
determining the number of dominant modes may succeed where peak counting 
falls. If such a method should reveal the presence of previously undetected 
modes, then careful curve fitting or some other modal testing technique, 
such as multiple-shaker tuning, might successfully separate the modes. 

The problem of determining the number of dominant modes was discussed 
extensively some years ago in connection with the number of shakers re- 
quired to separate modes in multiple-shaker modal testing. Traill-Nash [1] 
Introduced the "effective number of degrees of freedom" at a given frequency, 
which he defined as being the number of motion coordinates required to repre- 
sent with accuracy structural response at that frequency. He concluded that 
the number of shakers must equal or exceed the effective number of degrees 
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of freedom. Bishop and Gladwell [2 | suggested a relationship between Tralll- 
Nasfi' . < ffectlve number of degrees of freedom and the number of dominant 
modes; subsequently, Asher [3] Implicitly equated these two numbers. Me 
then stated. In effect, that the number of shakers required equals the 
number of dominant modes. This contention Is not generally true; the number 
of distinct generalized f^'ce distributions must equal the number of dom- 
inant modes, but there Is no necessary relationship between the number of 
generalized force distributions and the number of discrete forces. Nonethe- 
less, Asher made a significant contribution by proposing probably the first 
quantitative methods to determine the number of dominant modes by analysis of 
transfer function data. 

This paper describes the theoretical basis for such a quantitative 
method, reviews the methods discussed by Asher, proposes a new and 
more effective method, and illustrates these methods with the use of numer- 
ically simulated transfer function data. 

II. THEORETICAL BACKGROUND 

Consider a linear structure discretized to n degrees of freedom, the 
time-dependent responses of which are elements of the n x 1 column matrix : . 
(Notation Is listed at the end of the paper.) The governing matrix equation 
of motion is 

[m] x + [c] x ♦ [k] x ■ f (1) 

where [m], [c], and [k] are the n x n inertia, damping, and stiffness matrices, 
respectively, and f is the column matrix of time-dependent forcing. We 
specify that all forces vary sinusoidal ly at the same frequency, u, and that 
all have 0° or 180° phase, 

f * F cos wt 3 Re Jf e iu,t | (2a) 
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After starting transients have deca d away, response Is steady-state slnu 
soidal , 

x - Re } < e 1wt } (2b) 

where the elements of 'amplitude vector x are generally complex, reflecting 
phase differences between excitation and response. The linear, frequency- 
dependent relationship between excitation amplitude and complex response 
amplitude Is defined by the n x n transfer function matrix [HU)], 

X - [H( w )] F (3) 

The standard real modal analysis solution of equation (1) for [H(w)] 
begins with calculation of the real undamped natural frequencies w , r ■ 1, 

2, •••, n, and the associated real mode shape vectors * r « which are the 
columns of modal matrix [«♦,] (Mel rov Itch [4]). Subsequently, response co- 
ordinates are transformed Into normal coordinates which diagonalize the 
mass and stiffness matrices of equation (1); then the normal coordinates 
are calculated by matrix Inversion, and X Is calculated from the normal co- 
ordinate solution in the form 

X - [*] [SU)] F 

where [SU)] Is an n x n complex matrix. Hence, the transfer function matrix 
Is 

[HU)] - U] [SU;] 

and any column, say the jth, of [H(u>)] can be written as 

n 

HjU) ■ y S r jU) $r* J * l,2,***,n (4) 

r = 1 

Thus, each column of the transfer function matrix can be expressed as a sum- 
mation of the n linearly independent mode shape vectors. If we consider some 
subset p < n of degrees of freedom and define the corresponding p x 1 incomplete 
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Ith transfer function column as H. (j), then equation (4) gives 

p 

HjU) - ^S rj (u.) a* . j • 1, 2 n (5) 

rT 


where the degrees of freedom Include' In each p x 1 Incomplete mode shape vertoi 
{*, are the same as those Included In Hj . Although the summation In equation 
(5) extends over all modes, only p of the n vectors are Independent. 

If damping matrix [c] were to couple the undamped normal modes (l.e., If 
[*]* [c] [*] were not diagonal), then the use of undamped normal modes as out- 
lined above would be computationally Inefficient, and we would probably solve 
for the transfer function matrix with an appropriate complex modal analysis 
employing complex eigenvalues and node shape vectors. Nevertheless, the so- 
lution for each column of the transfer function matrix could still be expressed 
In the forms of equations (4) and (5), that Is. as linear summations of n mo te 
shape vectors, where In this case the $ r would be complex vectors. The Im- 
portant fact, expressed In the language of matrix theory. Is that each transfer 
function column H. Is an element In the n-dimeislonal vector space spanned 
by the n mode shape vectors, whether they be real or complex; similarly, each 
Incomplete column Hj Is an element in the p-dimenslonal vector space spanned 
by any p linearly Independent Incomplete mode shape vectors. 

A useful general characteristic of structural dynamic behavior is that 
very few of a structure's many vibration modes are sensitive to excitation at 
any given frequency. These few modes then dominate the response at that frequency. 
If there are q such dominant modes at frequency u, then the mathematical state- 
ments of their dominance, from equations (4) and (5), are 
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( 6 ) 


H j(w) S r j(ui) J " It 2# • • • |H 

q 

HJU) : V S rj U) ♦* . J ■ 1, 2 n (7) 

q 

whe^e £ denotes summation over only the q dominant modes. In equation (6). 
q 

the q n x 1 mode shape vectors associated with the dominant modes generate « 
q-dlmenslonal subspace of the original n-dlmenslonal vector space. The ap 
proximate equality In equation (6) means that each Hj(u>) column can, with 
small error, be considered an element of the subspace. In equation (7), the 
p x 1 Incomplete mode shape vectors associated with the dominant modes generate 
a q- dimensional subspace of the original p-dlmenslonal vector space, provided 
that p^q. Again, the apprr-. .mate equality means that each H* (u) column Is 
approximately an element of the subspace. 

Most current methods of modal testing are capable of measuring Incomplete 
transfer function column vectors over a frequency band of Interest. The jth 
column Hj represents physically the complex response amplitude of motion sensor?, 
at p stations on the test structure due to sinusoidal forcing excitation of 
unit amplitude at the jth station, which does not necessarily coincide with 
any of the motion sensor stations. If excitation Is applied successively 
to k different stations, then vectors H^j , j * 1, 2, ...k, can be measured. 
They are the columns of the p x k Incomplete transfer function matrix [H ]. 

(It is generally Impossible to measure the complete matrix [H], since a con- 
tinuous structure has an Indefinitely large number of degrees of freedom.) 

If the structure being tested responds linearly, then each H? column 
Is represented mathematically by equation (7), which, therefore, is the basis 
of the methods discussed in Section III below for determining the number of 
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dominant modes from experimental tr sfer function data. Each method estl 
mates the number q of dominant ve tors a* In equation (7), gl »n k ex 
perlmental p x 1 vew.ors HT . 

One practical requirement for the correct use of equation (7) In the 
present context Is Immediately evident: since p must exceed q, as dlscusvl 
above, the test engineer must guess an upper bound o for the number of 
dominant modes likely to be encountered, and then he must Install more than 
q distinct motion sensors. This requirement does not present a significant 

M X 

practical obstacle, since q for most structures should be on the order of 

nio x 

ten or less. It Is assumed In the remainder of the paper that the number of 
motion sensor measurements available for analysis Is always greater than the 
number of dominant modes. 

In vehicle modal testing. It Is usually feasible to install a substantial 
number of motion sensors, but the number of exciters or excitation stations Is 
often much smaller due to practical limitations. Hence, we assume that k p 

In most of what follows. 

III. QUANTITATIVE METHODS FOR DETERMINING THE NUMBER OF DOMINANT MODES 
III .1 TRANSFER FUNCTION DETERMINANT METHOD 

This method Involves analysis of square transfer function matrices, which 
are formed by the use of only k of the p available motion sensor measurements 
Thus, [H*] is a k x k matrix. If, In the first case, the number of exciters 
Is less than or equal to the number of dominant modes, k < q, then accordlnq 
to equation (7), the k columns of [H*] generally will be linearly independent: 
hence [H* ] will be non-singular and its determinant will be non-zero. 
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det IH*] / 0 for k S q. But \ r th« number of exciters Is Increased until if 

just exceeds the number of dominant 'nodes, k ■ q ♦ 1, then the k columns 

will be approximately dependent, so that [H*] will be nearly singular and 
Its determinant will be close to zero, det [H*] • 0 for k > q. The strategy 

for application of this method, therefore. Is to add rows and columns of la 

to the transfer function matrix In unit steps until the value of k Is found 
for which det [H*] *0; then the conclusion Is that q ■ k - 1. 

This Is a very simple method to apply, but It has some deficiencies. 

First, transfer function determinants are complex, so one must assess the 
possibly non-monoton!c progression toward zero of a sequence of complex numbers.- 
Second, the restriction to square transfer function matrices Is undesirable 
because It prevents the use of all available motion sensor data. Both of 
these deficiencies are eliminated with little additional effort by the use 
of the * • determinant method described below, so there appears to be no 
reason to develop further or test the transfer function determinant method. 

We note that Ibrahim and Ml kul elk [5] employed a similar method, but with 
filtered transient response data, and found It quite satisfactory . 

III. 2 GRAM DETERMINANT (GRAMIAN) METHOD 

This method Involves analysis of a rectangular p x k transfer function 
matrix, [H*]. The Gram matrix of [H*] Is defined to be a matrix product, 

[G k ] • [fl*] t [H*] 

where the overbar indicates complex conjugacy. By this definition, the Gram 
matrix Is a k x k Hermltlan matrix. The Gram determinant, or Gramlan, Is 
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defined to be 

G (k) . .let [G k ] 

It can be proved that the Gramlan Is real and non-negative. 

The Gramlan of a transfer function matrix Is a quantitative measure of 

the degree of linear dependence of the column vectors, Hj, j * 1, 2, 

Specifically, the set of vectors Is linearly dependent If and only If the 

Gramlan Is zero (Hildebrand [6]). Moreover, It Is reasonable to expect 

that if the set of vectors Is almost but not precisely linearly dependent, 

then the Gramlan should De nearly zero. As Is discussed In Section III.l, 

any q or fewer columns of [H*] generally are linearly Independent and, hence, 

Ik) 

have non-zero Gremians, G' ' > 0 for k ■ 1, 2, ...,q. But any set of more 

than ^ columns will be approximately dependent and have very small or zero 
(kl 

Gramlans, G' ' = 0 for k >_ q + 1 . Therefore, the basic strategy for appli- 
cation of this method Is to add columns of data to the transfer function 

(k) 

matrix in unit steps until the value of k Is found for which G x ' * Os then, 
the conclusion is that q = k - 1 . 

It Is necessary in applying this method to separate the change in Gram 

determinant value due to change In degree of linear dependence from that due 

simply to change in determinant size. If, for example, all Gram matrix elements 

_2 

are numerically o f order 10 , then, without change in degree of linear de- 

pendence, G^ ^ will be of order 1C G^ of order 10’^, of order 10 r \ 
and so forth. This characteristic of determinants will obviously mask the 
Gramian test for linear dependence unless Gram matrix elements are of order 1 
In applications of the Gramian method, the authors have attempted to minimize 
the masking by normalizing each column of [H*] so that the corresponding diaaonal 
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element of the Grom matrix Is 1, 1.e. t (H^ • M j ) " 1, J ■ 1,2,. ...k, where 
H* here denotes the normalized col urn rather than the original column In 
physical units. The numerical results of Section IV. 3 suqgest that this 
ad hoc measure Is effective In filtering out Gramlan variation with Gram 
matrix size. 

Asher [3] described and discussed both the transfer function determinant 
in l the Gram determinant methods. He recognized that In applying either method, 
one might find It difficult to decide how small a generally non-zero determinant 
value must be In order to Indicate linear dependence correctly. The numerical 
results of Section IV. 3 confirm that the absence of a definite smallness 
criterion is Indeed a weakness of the Gramlan method. Even though all Gramlan 
values are referenced to G^ ■ 1 by virtue of the normal Izatlon procedure de- 
scribed above, examples for different situations show Gramlan values cf 

-1 -2 -3 

orders 10 ,10 , and even 10 . 

Another deficiency of the Gramlan method is that it can correctly Indicate 
linear dependence, yet underestimate the number of dominant modes. To under 
stand how this can happen, consider a simple example in which there are two 
dominant inodes . The analysis of three given, distinct transfer function vec- 
tors, denoted Vj, v.,, anJ v^, then should indicate two modes. The set of 
three vectors is linearly dependent, but Suppose also that and v, are 
independent and V| and v ? are dependent. If the transfer function matrix is 
defined as [H*] » (.Vj, v,, v^], then G^ ■ 1, G^ > 0, and cf^ r 0, leading 
to the correct conclusion that q 3 2. If, on the other hand, we define [H*] 

LV| , Vj, v.J, then 6^ » 1 and G^ ■ G^* 0, leading to the incorrect con- 
clusion that q * 1 . An instance of this particular case occurring in a 
realistic situation is presented in Section IV. 3. It is clear that the 
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ordet ing of vectors In the p x k transfer function matrix affects all Gran » .w* 
values except ■ 1 and G^, wh r.h Is Invariant with column and row o» ,,, ''1'" 
The Gramlaii method then has some serious weaknesses. Perhaps for this 
reason. It apparently has not been employed widely. The authors have lo- 
cated only one published application, that by Klosterman [7], The vector- 
fit method to be discussed next is, to a considerable extent, free of the 
weaknesses of the Gramlan method. 

1 1 1. 3 VECTOR-FIT METHOD 

A concept analogous to the number of dominant modes of a vibrating struc 
ture is that of a "best approximating subspace". Cliff [8] discussed this 
concept in the context of control theory. Given a set of k p-dimensional 
vectors, one can calculate the particular m-dimensional basis (m < k) which, 
among all possible m-dimensional bases, does the best job of approximately 
spanning the set of k vectors, with error minimized In the least-squares sense. 
In other words, the k vectors are "fit" to the best approximating m-dimensional 
subspace. The method developed here to solve for the number of dominant modes 
follows Cliff's general approach; hence, it is referred to as the vector-fit 
method. 

Given the p x k transfer function matrix [H*] for frequency u, the general 
stepwise procedure for application of the vector-fit method Is as follows: 

1. A particular p-dimensional complex unit vector u^ is calculated from 
LH*]. Among all possible unit vectors, Uj alone has the property of 
producing the best, in a sense to be defined, set of one-tern approxi- 
mations to the transfer function columns. This set of appro<imations 
takes the form 

Hj » Cji u r j = 1 ,2, - • • , k 
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where the t^. generally are complex constants. Next, the real scalar 
error associated with this set of approximations Is calculated. 

2. A second unit vector u? Is calculated. It Is orthogonal to U|. Among 
all possible unit vectors orthogonal to u^ , u, alone has the property 
of producing, in conjunction with u, , the best set of two-term approxi- 
mations to the transfer function columns. This set takes the form 

2 

5 } Up j ■ 1,2,**», k 

1 » 1 

12) 

Next, error E' ' associated with this set of approximations Is calculated. 

3. The mth unit vector u Is calculated. It Is orthogonal to all other u., 

1 * 1,2,***, m-1. Vector u (]) has the property of producing, in conjunc- 
tion with iip u., ,••*, u jn _p the best set of m-term approximations to the 
transfer function columns. This set takes the form 

m 

HJ * C j^ Up j ■ 1,2,*“, k 

1 « 1 

Next error E^ associated with this set of approximations is calculated. 

Each step introduces a refinement of the approximation, so the error diminishes 
In each step, E^ 1 "^ ^ E^ m If, after m > 1 steps of this procedure, we find 
that E^'^i 0 relative to E^, then we may reasonably conclude that the set 
ol transfer function vectors is spanned approximately, with very small error, 
by an m-dimensional basis. According to equation (7), then, there are m 
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dominant modes at frequency u, i.e. q ■ m. It is quite unlikely that the '.nit »< u 
will he identical to the mode shape vectors ^ of equation (7), since the 
Uj are orthogonal by definition while the need not and generally will 
not be orthogonal. Nonetheless, it is certain that both sets u. and ** span 
tlie same q-dimensional vector space. 

If p < k and the procedure Is carried through p steps, then E ^ 0 

since the p independent u^ exactly span p-dimensional space. If, on the othn- 
hand, k < p and the procedure is carried through k steps, then 0 

since the k independent u. exactly span the subspace defined by the k trans- 
fer function column vectors. (This case is simply orthogonal Izatlon of a 
vector set, similar to the Gram-Schmidt procedure.) Hence, in order for the 
vector-fit method to produce a correct evaluation of the number of dominant 
modes, it is necessary that p and k each must exceed q. That is, the test 
engineer must provide both more motion sensors and more excitation stations 
than the maximum number of dominant modes likely to be encountered. In appl i 
cation, the method itself indicates if too few motion sensors or excitation 
stations have been used, and this is illustrated in Section IV. 3. 

It is worthy of note that each ot the three methods described for i den t i 
fication of the presence of q dominant modes requires a minimum of q + 1 motion 
sensors and q *■ 1 excitation stations. 

The theory associated with the calculations discussed above is developed 
next. Let a basis for p-dimensional space of complex vectors consist of p 
orthogonal unit vectors u^ , u. ; , •••, u () , which are unknown at this point. Any 
transfer function column vector can be expressed as the summation 


/ P \ 

HE* E Mv j 

'i = 1 i = m + 1 / 









The basis vectors are orthogonal In t ie Hermit Ian sense, so 

(fj- »j)-(E * E )(“r «J)(?3‘ “() 

\l ■ 1 1 ■ m ♦ 1 / 

Hence, we may define as follows the squared scalar error associated with 
approximation of H* as a series sum In only the first m basis vectors, 

M 2 ■ 2 Si) 

' ' 1 ■ m ♦ 1 

m 

■(?}• -j) ’ S (-r *i) 

i - i 

The appropriate total or global error, defined in the least-squares sense 


over all Hj, Is 


(m) . 


t (•'?) 

j 3 1 ' t 


2 - 

i 2 


r> in -j 

■ !J) - E E (sr si) 2 

j 3 1 j ■ 1 1 ■ 1 

lor given m, we wish \.o determine the basis vectors u^ so as to minimize E^ 1 "^ 
by maximizing the double-summation term. With a few steps of matrix algebra, 
that term is cast into the more useful form 


E E (V SjHSj- St) ■ E Si s, 

j s 1 1 3 1 1 3 1 


r 


where [A] ■ [H*] [A*] 1 Is a p x p Hcrmltlan matrix with p real non-negative 
eigenvalues x^ >_ _> . . . > x >0 and p corresponding complex, mutually 

orthogonal unit eigenvectors • i> 0 , .... . Thus, for given m, we wish to 

| k P 


minimize 

E (m) 



1/2 


yf [A] Uj 


Consider first m * 1 .■ The maximum value assumed by the Hermitlan quadratic 
form u ^ [A J u^ Is equal to the largest eigenvalue x^ , and this maximum results 
if Uj ■ ^ (Franklin [9]). Hence, the minimum associated error Is 


E o) 






1/2 


Next, for m = 2 we wish to minimize 


E 


( 2 ) = 





1/2 


For all possible u^ orthogonal to u^ * ^ , the maximum value assumed by quad 
ratlc form u^CA] u 2 is equal to the second eigenvalue X^. and this maximum 
results if u^ = ^ (Franklin [9]). Thus, 


e ( 2 ) 


k 



j=l 



1/2 


This reasoning can be extended easily to show that the minimum error for arbitrary 
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i r. ill 

E (m) . £<?;.«;). £*, 
j-i i«i 


where u^ ■ To simplify notation. It Is understood In equation (8) and In 

/ m \ 

the discussion and numerical examples to follow that E' ' denotes the minimum 

error, even though subscript min Is deleted. 

We can now reiterate and summarize the vector-fit method for determlnino 

the number q of dominant modes at a given frequency. From the p x k transfer 

function matrix [H*], the p x p Hermltlan matrix [A] ■ [H^HR*]* Is formed. It 

Is necessary that p > q + 1 and k > q + 1. Next, the real eigenvalues of 

[A],Xj > Ag > ... > Ap > 0, are calculated. Next, for m * 1, 2 p, the 

minimum error values are calculated from equation (8). The smallest Integer 

m for which E^ : 0 relative to E^ then Is equal to the number q of dominant 

modes. Note that the eigenvectors of [A], ^ ^ $ . need not be calculated. 

It can be proved easily that th? error values E ; are entirely Independent 

nf the numbering or ordering schemes used to identify motion sensors and exciters 

or to arrange elements in the transfer function matrix. In other words, any or 

ill of the rows of [H*] can be interchanged and/or any or all of the columns can 

(in) 

he interchanged without changing the values of E' . Hence, the vector-fit method 
has no weakness comparable to the dependence on vector ordering of the Granilan 
method. 


IV. NUMERICAL SIMULATION STUDY 

Two mathematical models have been used for numerical simulation of experi- 
mental application of the vector-fit and Gramian methods. The models were desiqned 
Hi have frequency bands of high modal density with prescribed numbers of dominant 
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modes. The basic objective of the numerical study then was to determine If the 
vector-fit and Gramlan methods are capable of correctly determining the number 
of dominant modes. Other objectives were to compare the two methods and to 
develop guidelines for applying the methods and Interpreting the results. 

To simplify calculations. It was specified that each model have hyste'etlr 
damping which does not couple the undamped normal modes. Hence, transfer function 
matrix elements were calculated exactly from the equation 


■u 


£ 

r«l 


'ir Mr 


pr 

2 

1 U) 

rr 

(1 - V } + % 


is the generalized mass of the rth normal mode, and ay and g r are, respectively, 
the natural frequency and hysteretlc damping of the rth mode. 


IV. 1 TEN-MODE MODEL 

The mode shapes of this model are those of the first ten out-of-plane 
vibration modes of the uniform stretched membrane shown on Figure 1, 


♦ir 3 A r sin 



where Integers ny and n^ are the numbers of half-wavelengths for the rth mode 
listed In Table 1, x^ and y. are coordinates of the motion sensor/exciter 
station.; shown on Figure 1, and constant for the rth mode is chosen so 
that the numerically largest mode shape element equals one. Other modal 
parameters, as listed In Table 1, are not those of the uniform membrane, but 
rather were selected to produce a mathematical model with four very closely 
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spaced modes in the vicinity of 100 on the frequency scale and six more modes 
outside the region of high modal density. Thus, q ■ 4 In a frequency band or 
roughly three units centered at 100, so the methods under study, if successful, 
should indicate accordingly. Since this model has limited physical significance, 
ill quantities are considered to be dimensionless. 


IV. 2 FIFTEEN-DEGR'.E-OF-FREEDOM MODEL 

This model, un ike the ten-node model, Is based entirely on a physical 
structure, the cantilevered, rectangular plane grid shown on Figures 2a and b. 

I he model was designed to have two modes with nearly identical natural frequencies, 
as shown In Table 2. To achieve such close modes, an optimization technique 
similar to that of Hallauer et al [10] was used. 

Each elastic member of the model is a steel bar having Young's modulus 
F ■ 200 GPa (29 x 10^ psi), shear modulus G ■ 82.7 GPa (12 x 10^ psi', and di- 
ameter of 6.35 mm (0.25 inch). (The fundamental units of pounds. Inches, and 
seconds were used In all calculations.) Each bar may twist about Its axis and 
bend out of the grid plane, and each is clamped at both ends by either the rigid 
support wall or a rigid cylindrical joint member of 0.051-m (2-1nch) diameter. 

The five node points of the model are the intersections of the bar centerlines 
in the grid plane. The fifteen degrees of freedom then consist of one out-of- 
rlane translation and two out-of-plane rotations of each node point. The trans- 
lations are identified in Figure 2a, and the rotations of node 1, for example, 
ire Identified in Figure 2b. The nodal lumped mass and moment of inertia (for 
both rotational dof) associated with each rigid joint are listed in Table 2. To 
account approximately for distributed inertia of the bars, finite element con- 

3 

sistent mass matrices were used, with bar density taken to be 7859 kg/m 







-4 2 4 

(7.3536 x 10 lb-sec /In ). Modal nysteretlc damping constants g r were 
specified to be 0.015 for modes 1 and 2, 0.018, 0.022, and 0.026 rc r ctlvely 
for modes 3, 4, and 5, 0.03 for modes 6-10, and 0.035 for modes 10-ia. 

IV. 3 RESULTS AND DISCUSSION 

Consider first the ten-mode model. Typical transfer functions are 
shown on Figures 3a and b, with the coincident (co) or real pans and the 
quadrature (quad) or Imaginary parts plotted separately. Qualitative In- 
spection of these and other transfer functions seems to Indicate the presence 
of only two dominant modes In the vicinity of frequency 100. But there are 
actually four dominant inodes In that region, as listed In Table 1. 

To begin numerical application of the vector-fit method, we examine 
first errors calculated from the complete 10 x 10 transfer function matrix. 
Such a large matrix will not generally be required, but examining this case 
first establishes a reference and provides guidelines for interpreting the 
error values and, more generally, for applying the method. Since transfer 
functions are usually plotted versus excitation frequency, it seems natural 
and. In fact, proves advantageous also to calculate and plot error values 

t m \ 

versus excitation frequency. Errors E v ‘ ni ■ 1,2, •••,10, for the complete 
transfer function matrix were calculated from equation (8) and are plotted 
on Figure 4. Note the relatively large variations with frequency and the 
peaking of E^, E^, and E^. Note also that E^ Is on the order of 10 
of the maximum value of E^, and that E^ varies only slightly with fre- 
quency. As stated in Section III. 3, the result that E^ << E^ indicates 
that there are four dominant modes. But It is also important that E^ f 0 

and that E^, E^, etc. vary slowly with frequency In comparison with E^ * 
(2) . (3) 

E , and E . From these observations, we conclude that four modes dominate 



but do not completely describe the response in this narrow frequency band, 
and that the contributions of distant 'nodes remain relatively constant In 
this band. Both conclusions are quite reasonable in view of the nature of 
the ten-mode model Thus, this reference case suggests that 1 examining 
graphs of error • us frequency, we should use not only the basic criterion 
^(q) << a lso additional criteria that E^ i 0, that E^. 

E^ * etc. should vary slowly with frequency In comparison with E^, 

E^, •••, E^ ’ and that plots versus frequency of E^, E^,**., E^' 1 
should exhibit peaks. An exceptional case for which these criteria might 
not apply Is that of a distant mode contributing significantly to response. 

In this case, we might expect to find a relatively large, slowly varying, 
and non-peaking E^ 1 "^ for m < q. 

Figures 5 a-d are plots of errors calcu.aed from responses at seven 

motion sensor locations due successively to excitation at three, four, five 

and six locations. This simulates a realistic testing approach In which 

very few exciters are applied Initially, and additional exciters are applied 

12 ) 

as required. In Figure 5a for three exciters, E 1 # Is not much smaller than 
E^, and E^ ■ 0, as required by the theory of Section III. 3. (The very 
small non-ze^ values on this :nd other computer-generated plots are due to 
round-off errors.) Clearly, more exciters are required to indicate correctly 
the number of dominant modes. The addition of a fourth exciter ( 1 . e . , another 
column In [H*]) leads to Figure 5b, which also Indicates the need for at 
least one more exciter. With five exciters, however. Figure 5c shows that 
E^ varies slowly with frequency and is non-zero yet much smaller than E^. 
We conclude, therefore, that q ■ 4. Figure 5d for six exc'ters substantiates 
the conclusion. 

Unfortunately, interpretation of the error plots Is not always as 
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unambiguous as It appears to bo with Figures 4 and 5. There Is a definite 
dependence on the motion sensor and exciter locations represented In the 
transfer function matrix. These points are Illustrated on the error plots 
of Figures 6 a,b, and c, which were calculated for different sets of notion 
sensor/exciter locations than those of Figures 4 and 5. Figure 6a for f 1 v»* 
notion sensors and five exciters can be Interpreted as Indicating four 
dominant nodes; but this Interpretation Is substantially weaker here than 
for Figures 4 and 5c, d since E ^ on Figure 6a Is only slightly peaked and 
Is very small relative to E^. Flyure 6b for six notion sensors and fix 
exciters Is ambiguous; one might Infer that It Indicates three, four, or 
five dominant nodes, with three being perhaps the most likely Interpreta- 
tion. Figure 6c, also with six notion sensor/exciter locations but one dif- 
ferent location than Figure 6b, permits a somewhat more certain Interpreta- 
tion of four dominant nodes. 

One can observe from the vertical axis scales on Figures 4-6 that 
all non-zero error values tend to Increase as columns are added to the 
transfer function matrix. This tendency appears not to hav*» any useful 
significance. 

Next, we examine the Gramlan method as applied to the ten-node model. 

Figure 7a Is a graph of Gratnlans plotted versus frequency in the region of 
high nodal density. Gramlans G^ - G^ 10 ^ were calculated from complete 10 x 1 
transfer function vectors, with the vectors applied In the order of the station 
numbering shown on Figure 1, l.e., 1, 2, .... 10. Recall that the identifi- 
cation criterion for this method is G^ + ^ 5 0. where G^ » 1 by virtue of 

(51 

transfer function vector normalization. So G' should be nearly zero in this 
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case. A tabulation of G x for Figure 7a (not Included) shows values In 

the ranges 10* 3 and 10* 4 . Figure 7b Is a graph of Gramlans G^ - G^ 10 ^ 

also calculated from 10 x 1 transfer function vectors, but with the vectors 

applied In Inverse or.der of the station numbering, l.e., .0, 9,...,1. Values 

of for Figure 7b are In the ranges 10* 1 and 10*^. The vast differences 

between the curves on Figure 7a and those on Figure 7b Illustrate the vectoi 

ordering dependence of the Gramlan method, and the different magnitudes of 
(5) 

G' on the two graphs Illustrate the Indefiniteness of the smallness criterion, 

G* q+1 * : 0. 

Finally for the ten-mode model, Figures 8a and b demonstrate the slgnU* 
cance on typical transfer function plots of the dominance of four modes. Ihes' 
graphs show exact coincident and quadrature values of j(w) and H 7 7 (w), aproxl 
mate values calculated from only the four dominant modes, and approximate values 
calculated from only modes four and six, two of the dominant modes. Whereas les 
four and six alone dominate 3 , all four of the close modes make comparable 
contributions to H ? y. 

Consider next the flfteen-dof model, which has the pair of close modes 
listed In Table 2. In testing of this model. It would be natural to Instrument 
and to provide forcing excitation at the five translation degrees of freedom 
shown on Figure 2a, so we will ana yze the 5 x 5 Incomplete transfer function 
matrix associated with those degrees of freedom. 

Tlgures 9a, b, and c are selected transfer functions In a narrow band about 
the close natural frequencies. Figure 9a for ^(w) has the character of n* st 
transfer function elements In this frequency hand, namely, It seems to In- 
dicate only a single. Isolated mode. Figure 9b for jqM suggests the 
presence of more than one mode, but the asymmetry of the coincident response 
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curve would probably be attributed lo a distant mode rather than to close modes. 
Of all the elements of the 5 x 5 [H*], only ^(w) shown on Figure 9c pro- 
vides definite evidence of the presence of two close modes. But both H 4 .10 

and Hj j-j have such small magnitudes relative to ^ that. In actual testing, 
they would probably be lost In noise or ignored. So It Is fair to say that 

qualitative examination of the transfer functions Indicates the presence of 

only a single mode at 149.2 rad/sec. 

Before examining error plots in the region of the two close modes. It is 
useful to have as a reference an example of error plots at and near a single 
Isolated node. The second mode of this model, with a natural frequency of 
67.5 rad/sec. Is quite distant from all other modes. Figure 10a Is the graph 
of error values around this mode for the 5x5 [H*]. For comparison. Figure 
10b is the graph of error values In a region of almost no modal activity be- 
tween the second and third modes. The error scales of both Figures 10a and b 
are quite small (relative to that of Figure 11a discussed below), and neither 
figure has any error peaks. The only significant difference in character be- 
tween the two figures is the numerical noise at and near the natural frequency 
in Figure 10a, due to accumulated round-off error in eigenvalue calculations. 

Figures 11a and b are graphs of error values in a band around the pair of 
close modes, the former for excitation at all five translation degrees of free- 
dom, and the latter for excitation at only three. On the basis of all criteria 
developed previously and in comparison with Figure 10a, these graphs indicate 
clearly and indisputably the presence of two dominant modes. 

Figures 12 a, b, and c are graphs of Gramians in the region of close modes. 
Figure 12a for five exciters applied in the order 1, 4, 7, 10. 13 seems to 
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indicate clearly that q 2. So also does Figure 12b for three exciters ap- 
plied In the order 13, 1, 10. However, Figure 12c for the same three exciters 
applied In the order 1, 10, 13 seems to Indicate Just as clearly, though In- 
correctly, that q ■ 1. 

In summary of the discussion of theory and the numerical simulation study, 
the vector-fit method Is distinctly superior to the Gramlan nx'thod for the pur 
pose of determining the number of dominant modes. The only advantage of the 
Gramlan method Is that it requires substantially less computation time. As 
has been demonstrated, the vector-fit method can produce error plots which ire 
difficult to interpret correctly. Hut results of the simulation study suggest 
that the use of a large number of motion sensors In calculation of the error 
values will reduce the likelihood of mistaken Interpretations. Theoretically, 
a minimum number of q+1 notion sensors are required; however, it would seem 
prudent and usually practical to estimate q and then to use several times that 
number of motion sensors In calculating error values. 

We note that the vector-fit method Is valid regardless of the type of motion 
sensor employed in testing, since the form of equation (7) remains the same fo> 
displacement, velocity, or acceleration transfer functions. 

V. CONCLUIUNG REMARKS 

The vector-fit method for determining the number of dominant vibration nod*’s 
from structural transfer functions has been derived theoretically, illusti ited 
with a numerical simulation study, and compared with other methods. The method 
works well with exact, noiseless, simulated data. However, its practical ap 
pllcabillty has not been evaluated. The logical next step is to test the method 
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with real data. If it should prove applicable, It would be a useful analysis 
tool for modal vesting. Regardless of the type of test 4 ng employed, the vector- 
fit method could provide an independent check on whether or not all significant 
modes had been de tec tied. 
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mth real error value 
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G'^ ■ det[G^] 
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Grain determinant, Gramlan 

hysteretic damping constant of rth mode 


number cf dof subjected to forcing excitation 

ith eigenvalue of [A] 

generalized mass of rth normal mode 

numbers of half-wavelengths for the ten-mode model 

total number of dof 

frequency of excitation 

natural frequency of rth normal mode 

number of dof Instrumented with motion sensors 

number of dominant modes 

element of [5] 
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TABLE 1 




Parameters of 

ten-mode model 



r 

1 2 

3 4 

5 

6 

7 8 

9 10 

'Op 

67.6 82.0 

91.2 98.9 

99.4 

100.1 

101.4 110.7 

117.6 132.3 

'V 

0.015 0.017 

0.017 0.021 

0.019 

0.023 

0.024 0.026 

0.027 n. 

'V 

0.02 0.009 

0.01 0.012 

0.011 

0.01 

0.009 0.011 

0.01 ’ 1 

r ,n r 

tt — f.r 

1.2 3.1 

2.2 

3.2 

1,3 4,1 

2.3 47 F 



TABLE 2 





Inertia parameters and close modes 

of fifteen-dof model 



Inertia parameters 



Close modes 



Mass Moment of inertia 


r 

3 

4 

Node 

kg 

10" 3 kg-m 2 

w 

r (r/s) 

149.18 

149.2? 

r 

(575TT4 

0.5570 


*i7 — 

1.000 

1 .000 

2 

1 .4994 

1.7199 


*4,r 

0.031 

-0.274 

3 

2.0390 

4.1627 


^7,r 

0.446 

-0.444 

4 

0.4756 

0.1154 


*10, r 

-0.421 

-0.236 

5 

2.0424 

4.1833 


*13, r 

-0.151 

0.216 
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Figs . 6 

Error graphs for the tui-morie model 
(a) motion sensors/exciters 2, 3, 4 
7 . 13; (b) motion sensors/oxci t«rs ? 
3, 4, 7, n, 10; (c) motion sensors/ 
exciters 2, 3, 4 , 6. 8 10 
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Figs. 8 Comparisons of exact and approximate transfer functions 

for the ten-mode model, with solid lines for exact values, 
long dashes for modes 4, 5, 6, and 7, and short dashes for 
modes 4 and 6: (a) 3 (u>); (b) H, 7 (w) 
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Figs. 12 

Gramlan graphs for the flfteen-dof 
model with motion sensors 1,4, 7. 
10, 13: (a) five exciters applied 

in the order 1, 4, 7, 10, 13; (b) 
three exciters applied In the order 
13, 1, 10; (c) three exciters applli 
In the order 1 , 10, 13 
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